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Abstract 

A simple analytical tool based on stochastic ordering is developed to compare the distributions of 
carrier-to-interference ratio at the mobile station of two cellular systems where the base stations are 
distributed randomly according to certain non-homogeneous Poisson point processes. The comparison is 
conveniently done by studying only the base station densities without having to solve for the distributions 
of the carrier-to-interference ratio, that are often hard to obtain. 

Index Terms: Carrier-to-interference Ratio, Co-channel Interference, Fading channels, Stochastic 
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I. Introduction 

A Poisson point process has been adopted in the literature for the locations of nodes in the 
study of shotgun cellular networks, ad-hoc networks, and other uncoordinated and decentralized 
' communication networks JTJ, [El and references therein]. An underlying assumption in all the 
previous work is that the density of transmitters, referred to as base station (BS) throughout 
this paper, is constant, i.e. the Poisson point process is homogeneous. Such a model does not 
sufficiently represent reality. 

In B2l and here, we have modeled BS arrangement by non-homogeneous Poisson point pro- 
cesses in M 1 , I = 1,2, and 3. We aim to characterize carrier-to-interference ratio (j) at a given 
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mobile station (MS). In Q, we have derived semi-analytical expressions for the tail probability 
of j, denoted as P({y >?/}■), by deriving the characteristic function of the reciprocal of j. 
Moreover, the j characterization holds for a wide range of scenarios of interest such as arbitrary 
distributions for fading and random BS transmission powers, arbitrary path-loss models and 
arbitrary locations for the MS. In spite of such a general result, it is still not convenient to 
compare the y distributions of two different networks. 

Is it possible to qualitatively compare two j distributions by only examining the BS densities 
without having to obtain the j distributions? This paper answers the question affirmatively 
for certain BS densities by developing a stochastic ordering based tool. Concepts of stochastic 
ordering have been applied to scenarios of interest in wireless communications in [3J. To the best 
of our knowledge, this is the first work that uses stochastic ordering to understand large scale 
random wireless networks. The main result of this paper is Theorem Q] for which Section Hill 
develops the necessary tools. The utility of this result is explored in Section [IVl by considering 
several scenarios of interest in modeling the wireless network. 

II. System Model 

The Shotgun Cellular System (SCS) is a model for the cellular system in which the BSs 
are distributed in a /-dimensional plane (I— D, typically I = 1,2, and 3) according to a non- 
homogeneous Poisson point process in M. 1 . The intensity function of the Poisson point process 
is called the BS density function. 

Without loss of generality, we restrict our attention to 1-D SCSs, because for the y analysis, 
the l-D SCSs can be reduced to an equivalent 1-D SCS [2, Lemma 2] with a BS density function 
A(r), where r > is the distance of the BS from a mobile-station (MS) located at the origin. 
For example, a homogeneous l-D SCS with density Ao (> 0) is equivalent to a 1-D SCS with 
density function A (r) = Aofyr' -1 , r > 0, b\ = 2, b 2 = 2tc, and 6 3 = Air Corollary 1]. 

The BSs are assumed to have independent and identically distributed (i.i.d.) random transmis- 
sion powers Ki's and shadow fadings ^//s across BSs. The deterministic path-loss is R~ £ , e > 0. 
We assume an interference limited system and omit background noise. We focus on the signal 
quality of a MS at the origin. The MS chooses to communicate with the BS that corresponds 
to the strongest received signal power, referred to as the "serving BS". All other BSs are the 
"interfering BSs". The signal quality at the MS is measured by j = ^, S ^^ S R _ £ , where S 
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indexes the serving BS and % indexes the interfering BSs. Further, R$ < -Ri < R% < • • • are 
ordered BS locations. 

III. The Stochastic Ordering of ^ 

In this section, we present the theoretical tools that are used to compare j tail probability by 
comparing the equivalent 1-D BS densities A(r). Since the effect of i.i.d. shadow fading factors 
and i.i.d. transmission powers can be captured by modifying the BS density as shown in Section 
|IV-D[ they are assumed to be 1 for all BSs. The generalization to arbitrary path loss model is 
given in [2, Section VI], which is also equivalent to modifying A(r). As a result, T = OQ 1 e . 

Definition 1. Let X and Y be two random variables such that P ({X > x}) < P ({Y > x}) , V x E 
(—00, 00), then X is smaller than Y in the usual stochastic order and this is denoted by X < st Y. 
Further, X = st Y means P ({X > x}) = P ({Y > x}) , V x G (-00, 00) . [4, p. 3] 

If X and Y are the j at the MS in two different SCSs, X < st Y implies that the MS in 
the SCS corresponding to Y is more likely to achieve better signal quality than in the SCS 
corresponding to X. 

Let {Rk}kLi represent the set of distances of BSs from the MS (indexed in the ascending 
order of the distance), D k+1 = Rk+i — Rk be the distance between two adjacent BSs, and 
fD k+1 \R k (d\ r; A(s)) be the probability density function (p.d.f.) of D k+1 conditioned on R k = r, 
as a function of the BS density A(s). 

Lemma 1. 



f Dk+llRk (d\r;X(s)) = e~K x ^ ds X(r + d), and (1) 
f aDk+1 \aR k (d' \r'; AGO) = f Dk+1 \ Rk (d' r'; -A (-) ) . (2) 



a \clj 

Proof: Equation ([!]) follows from the properties of Poisson processes [EJ , [I6J . Equation © is 
proved by f aDk+llaRk (d>\r>; X(s)) ( = } \f Dk+1 \ Rk ( d - 1-; X(s)) i ^A (^) exp (- X(s)ds^ 

= ^A exp ^— J^, +d ^X (^-) ds'^j , where (a) is obtained by a variable change; (b) follows 

from (OQ); and (c) is obtained by a variable change and gives ©. ■ 
Lemma Q] means that scaling D k +i and R k by a is equivalent to scaling the BS density as 
-A(-). The significance of Lemma \T\ in the context of j is as follows. 
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Corollary 1. The distribution of ~ at the MS in the 1-D SCS with BS density function A(r) is 
the same as that in 1-D SCSs with BS density functions ^A(^), V a > 0, i.e., yL n = s t y i\rzy 

Proof: Let {-Rfc}^ correspond to the 1-D SCS with BS density function A(r). Then, since 
the ordered BS locations R^s are determined by inter-BS distances, it follows from Lemma 

in that o\ = ( Q y;v. K: 

/ I A(r) E fc = 2 ( a - R fe) 



, where it^'s corresponding to -A( L ) 

a > a ' 



A(r) ESLa(^j 

have the same distribution as ai4's with A(r). 

As a result, {^A (^) , r > 0, a > 0} forms a parametric family of BS density functions 
such that the 1-D SCSs corresponding to them have the same j at the MS. In other words, 
appropriately scaling the BS density function will not change the p.d.f. of j. Moreover, the 
following special case is a direct corollary of the above result. 

Corollary 2. In a homogeneous l-D SCS, j is not a function of the BS density. 

Proof: Firstly, recall that the j at the MS in a homogeneous l—D SCS with BS density 
Ao is the same as that in a l—D SCS with a BS density function A (r) = Aofyr' -1 . Next, from 
Corollary Q] the distribution of ~ in this SCS is the same as that in a 1-D SCS with the BS 
density function ^A (-) = A a^r' _1 , a = a~ l , a > 0. Thus, distributions of y corresponding 
to aA and A are the same. ■ 
This was also observed in [O, where we showed that the expression for the characteristic 
function of (j) 1 did not involve A . Corollary [2] provides a simpler and more fundamental 
proof. Next, we define a notation used in Theorem [TJ 

Definition 2. For BS density function A(r), the cumulative BS density function is defined as 
/i(r) = J Q r \(s)ds, and its inverse function is define as [i~ l (q) = supjr : fi(r) < q}. 

Since A(r) > 0, /i (r) is a monotonically increasing function of r. In general, the inverse 
function is not injective since A(r) can be zero in arbitrary intervals of r e [0, oo) . The above 
definition makes it injective. For certain BS densities, it is possible to compare two j s by 
comparing the densities without solving for the distributions. It is facilitated by Theorem Q] 

Theorem 1. Let {Ai (r) , /i X (r) , /i^f 1 (g)} and {A 2 (r) , /i 2 (r) , fi^ 1 (q)} be the BS density func- 
tions, cumulative BS density functions and their inverse functions for two 1-D SCSs, respectively. 
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The ^ at the MS follows the stochastic order ^ L , , < st ~ I , , if for each q > and a = 

I I lAi(r) — / lA 2 (r)' * Ml ^g)' 

See Appendix [A] for the proof. Applications of the above theorem are in the next section. 

IV. Applications of the ^ stochastic ordering 

A. Comparison of Homogeneous l—D SCSs (I = 1,2, and 3) 

Here, we show that the signal quality degrades as the dimension / of the homogeneous l-D 
SCS increases, for which we need the following corollaries. 

Corollary 3. For each q > and a = if AA(r) = ^Ax(^) — A 2 (r) is a non-decreasing 

function for all r > 0, then tL , , <st tI, , 

— ' / lAi(r) — Sl / lA 2 (r) 

^{q) H^il) ^2 1 (<l) 

Proof: Note that f ^Ai (-) ds = q = f X 2 (s)ds. Hence, f AX(s)ds = 0. Suppose 
o a a 

AA(/i^ 1 (q 1 )) < 0, then AA(r) < 0, r G [0, /i^ 1 (<?)]> since AA(r) is non-decreasing. This is a 

contradiction. Thus, AA(/U^" 1 (g)) > 0. Using Theorem [U the corollary is proved. ■ 

Corollary 4. For a homogeneous l-D SCS with BS density Ao and its equivalent l-D BS density 
function A; (r) = Ao^r' -1 , r > 0, multiplying A; (r) with a non-increasing function f3{r) > 
improves the j, i.e., j| A;( - r ) <st jI^u^)- The inequality reverses if /3(r) is non-decreasing. 

Proof: If /3(r) is non-increasing, for any a > 0, the density difference AA(r) = ^A;(^) — 
f3(r)\t(r) = X^bi - I3{r)) r 1 - 1 is non-decreasing. By Corollary [31 j\ Xi{r) <st j\ mMr) 
holds. If j3(r) is non-decreasing, the same proof applies with AA(r) = /3(r)A;(r) — aA^(-). ■ 

(a) (b) 

Hence, j| Ai(r) >st f| A2(r) >st f| A3(r) b Y Pegging I = 1,2, 3 in A, (r) , respectively; (a) holds 
because A 2 (r) = /3 (r) Ai (r), where /3 (r) = is a non-decreasing function; and similarly (6) 
also holds. Thus, the comparison between — s is done without finding their distributions. 

B. A Qualitative Comparison between Two l-D SCS 

Consider a homogeneous l-D SCS with a BS density function X 1 (r) = A, r > 0, and 

{a < r < p 
, where a > /3. Such 
/3 r > p 

A 2 (r) might describe, for example, a highway passing through a region of greater population 
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(BS density of a) and then a region of smaller population (BS density of (3). Usually, such a 
scenario is approximated by a constant BS density through out the highway, which is represented 
by A x (r) . If A = a, it is easy to guess that j| Ai( - r ) <st t|w,.v ^ ut ^ ^ > a ' ^ * s not c l ear 
which SCS has better j. Theorem Q] shows that the SCS with density A 2 (r) has better j than 
Ai (r) , and the result holds irrespective of the specific values of a, (3, and A. To apply Theorem 

{I ' q<ap 

. Further, a(q) = 

^ , g < ap, I a > a = A 2 (r) , ^ < r < p 

. As a result, for q < ap, -Ai(-) = < , 

1+(/ ^7" )p/g ,g>«p " " [a>(3 = \ 2 (r) ,r>p 

and for q > ap, p 2 ^g) > P, ^i(£) = i+tf- a ) P / q > = r > p. Thus, applying Theorem 

EtLm ^ f| A2(r) - Similarly, if a < 0, we can show f | Mr) > st f ^ . 

C. Comparison of Path-loss Models 

Here, we compare the j at the MS in two homogeneous l—D SCSs with a BS density A 
and with different path-loss models, j^jyj, and j^yj, V r > 0. The proofs for Corollary \5\ and 
[6] involve the following common steps. Firstly, using Lemma 2], reduce the homogeneous 
/-D SCS to the equivalent 1-D SCS with a BS density function \(r) = Xokr 1-1 , i = 1, 2. 
Next, using [2, Theorem 5], the resultant 1-D SCSs can be reduced to an equivalent 1-D SCS 
with path-loss exponent e — 1, i.e., = ~, and with BS density functions Aj (r) = ^j^rj^y, 
z = 1,2 where /i^ 1 (■) is the inverse function of hi (■) . Finally, we have two 1-D SCSs with the 
same path-loss model = ~ and BS density functions A* (r) , r > 0, z = 1, 2. We show that 
Theorem Q] applies and derive the result. 

Corollary 5. In a homogeneous l—D SCS, if the path-loss follows a power-law parametrized 
by a path-loss exponent, e, the y at the MS improves as the path-loss exponent increases. In 
other words, if hi (r) = r £ \ i = 1,2, such that ei > e 2 > then (j) > st (j) , where (j). 
corresponds to the path-loss model 

Proof: At the end of Step 2, the equivalent 1-D SCSs with a path-loss model - have the BS 
density functions \ (r) = ^-r" 1 . Further, A 2 (r) = /3 (r) Ai (r) , where /? (r) = g-r^a - ", r > 
is a non-decreasing function. Hence, Corollary |4] applies and j\^r r \ >st jlwr) ' B 
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Hence, a simple proof that does not require solving the distribution of j gives the expected 
result that a channel with a greater path-loss exponent has a better j. The following corollary 
establishes a similar result between two popularly used path-loss models Q. 

Corollary 6. In a homogeneous l—D SCS with a BS density Ao, the received signal of a MS 
located at the origin satisfies < st (y) 2 , where corresponds to path-loss with 

{r £l , r ^ 1 
r £2 , r > 1 

where £ 2 > £i > /. The opposite conclusion holds when E\ > e 2 > L 

Proof: At the end Step 2, Ai (r) = ^-r~ 1 , r > 0, and A2 (r) satisfies the equation 

{1 , r < 1 

; ; . Since e 2 > £1 > I, (3 (r) is a non 

£i r ^2~~n , r > 1 

increasing function. As a result, Corollary |4] holds and hence y|w r -j <st 7\\ 2 r r y Thus, the 
system with the path-loss model has a better signal quality compared to that of j^^- Now, 
when £1 > £2 > /, p (r) is a non decreasing function and ylw^ >st T I A 2 (r) ' ' 

D. Shadow Fading and Random Transmission Powers 

In all the results until now, the shadow fading factors and the transmission powers for all 
the BSs were constant. Here, we generalize to the case when they are random variables, i.i.d. 
across BSs. The transmission power and the shadow fading factor of the same BS could be 
dependent. The following result reduces the 1-D SCS with random shadow fading factors and 
random transmission powers to a 1-D SCS where both are 1. 

Theorem 2. For the j analysis in a 1-D SCS with BS density function A (r) , if the random 

shadow fading factors {^i}^ and transmission powers {K i ] ( ^ =l are i.i.d. across all the BSs, the 

r 1 / 1 

SCS is equivalent to another 1-D SCS with BS density function A (r) = E$,,k (^K) e A ( r (^K) e 
where E is the expectation operator w.r.t. \l/ and K, which has the same distribution as and 
Ki, V i, respectively. This holds as long as the expectation converges. 

Proof: For the random shadow fading and transmission powers case, the j defined in 

Section UJ can be written as — = — ^— y e , where the index j corresponds to the BS 

E (RiiK^i)--) 
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with the strongest received signal power at the MS. The above expression for the j corresponds 
to a 1-D SCS with distances from the MS given by (K^i) £ j , unity shadow fading 
factors and unity transmission powers at each BS. Now, B2l Theorem 4] applies, and the new 
1-D SCS has a BS density function A (r) . ■ 
Thus, if there are random shadow fading factors and/or transmission powers at each BS, 
one can first apply Theorem [2] to obtain the equivalent 1-D SCSs with constant shadowing 
and transmission powers, and then apply Theorem Q] for the comparison of j s. The following 
corollary shows a scenario where ~ distribution is unaffected by the distributions of random 
shadow fading factors and transmission powers. 

Corollary 7. In a homogeneous I— D SCS with a BS density A , the ~ distribution at the MS does 
not depend on the distributions of the random shadow fading factors {$>i}°l 1 and transmission 
powers {K i \^t l , if they are i.i.d. across BSs and E# K (^/K)~ < 00. 

Proof: Firstly, recall that the homogeneous /— D SCS is equivalent to a 1-D SCS with a BS 
density function A (r) = X kr l -\ r > 0. We have (f)^ ( = } st f | X(r) ( = } st j\± x f±\ =st j| A(r) , 
where (a) is obtained by applying Theorem [2] to the 1-D SCS with the BS density function 
A (r) , to obtain the equivalent 1-D SCS with constant shadow fading factors and transmission 
powers and with a BS density A (r) = E*,x (9K)< A (r (VK)^ = E^ (VK)* A (r); (6) 

is obtained by rewriting A (r) as ^A (^) where a = ^E^^ (tyK)* j ' ; (c) is obtained by 
applying Corollary [U Thus, the shadow fading and random transmission powers have no effect 
on the j distribution. ■ 
This result was already proved in [2, Remark 4(a)]. But here, we have shown an elegant alter- 
native proof that is based only on the concepts of stochastic ordering. Finally, as a consequence 
of Corollary H the results in Section ITV-CI also hold for cases with random shadow fading factors 
and transmission powers that are i.i.d. across BSs. 

V. Conclusions 

This paper is an extension to our previous work in characterizing the j of a SCS in ||2). 
The study of the j at the MS in a cellular system with BSs distributed according to a non- 
homogeneous Poisson process is difficult because the distribution of the j is not in closed form 
and it is difficult to form an intuition about such networks. As a result, most of the j analysis 

October 17, 2011 DRAFT 



9 



are restricted to the homogeneous I — D SCSs. Here, we have developed a stochastic ordering 
based tool to analyze the y at the MS in such non-homogeneous Poisson processes. Due to 
Theorem [H for certain BS densities, we show that, by just comparing the BS density functions 
of the SCSs, we can make strong inferences such as, a MS in a given SCS achieves a y that is 
at least as good as that achieved in another SCS without having to solve for the y distributions. 
Moreover, as a consequence of Theorem [TJ elegant proofs are derived to show that (1) a MS 
in a homogeneous l-D SCS sees decreasing signal quality as dimension / increases; (2) the y 
at the MS improves as the path-loss exponent of the channel increases; and (3) as far as y is 
concerned, a SCS with random shadow fading factors and random transmission powers, which 
are i.i.d. across BSs, is equivalent to a SCS with constant shadow fading factors and transmission 
powers and with a modified BS density. 



Appendix 

A. Proof of Theorem [7] 

Consider the l-D SCS specified by the set {A (r) , /i (r) , n~ l (q)}, as in Definition [2l The 
following remark relates y to the cumulative BS density. 

If Ri denotes the distance between the serving BS and MS in the l-D SCS, 



if ( { j > y 



(a) 



ri=0 



( J > y 



Ri=r) f Rl (r)dr 



(b) 



q=0 



Q = q) fcj(q)dq, 



where Q = n(R\), and Q is an exponential random variable with mean 1. 

Equation (a) is obtained by conditioning w.r.t. R 1 . Equation (b) is obtained by express- 
ing (a) in terms of Q, where the p.d.f. of R± at R\ = (g) is f Rl (r)dr 



\r=li- 1 {q) 



-/o r ^)^A(r)(ir 



r=n 1 {q) 



e q dq = fQ(q)dq, which does not depend on A(r) 



To show that the BS density Ai(r) gives a worse y than A 2 (r) does, one needs to show that 



c 



< 21 



Tic w i s -m ^st 7i D -i/- -s w i ^ -V f° r a U <? > 0' where the condition of 

/ \R 1 =H 1 (q), Ai(r), r>/Ltj (q) — I \Ri=[i 2 (q), A 2 (r), r>^ 2 (q) * 

the domain of the BS density is because the locations of interfering BSs only depend on the BS 



density in that domain. Next, define a 



^2 *(g) 



. By Corollary [Q 



C 
1 



Rl^^iq), Ai(r), r>^ x ( g ) 



Eoo p_ 
k=2 n k 



R2>^ 1 (q), Ai(r), r>^-\q) 



H 2 >Mi (9). Ai(r), r>^ x (q) 



Y,T=2 ( R 'k) 
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where ii^'s are the ordered BS locations of the SCS with BS density ^A (0. The equation 
means that the conditional j of the SCS with a BS density Ai(r) is equivalent to an SCS with 
BS density ^Ai(^) with the same location of the serving BS as the SCS with BS density A 2 (r). 

With the locations of the serving BSs equal and fixed, j is a decreasing function of the 
interference. Theorem l.A.3.(a) of [4] says that decreasing functions reverse stochastic order. 
Therefore, one only needs to show that the interferences satisfy 

oo oo 
k=2 k=2 

As shown in [0 Appendix B], the total interference power can be expressed as YlT=2 R 'k £ = 
lim rB _ s . 0O limAr^oo J2f =2 Xi, where Xi is a Bernoulli random variable defined by 

F({X i = 0\R 1 =r 1 })=p il F({X i = r; £ + o(Ar)\R 1 =r 1 }) = l-p h 

Pi = A(rj)Ar + o(Ar), = r\ + (i — l)Ar, Ar = rB ^ ri , and r\ = ^ Now, since the 
condition ^Ai(^) > A 2 (r) holds for all r > /i 2 1 (0')' we have Xi|i Al (r^ > st Xj| A2 ( r ), Vi > 2. 
Since summation preserves the stochastic order [|4l Theorem l.A.3.(b)], © is proved, completing 
the proof of Theorem [Q 
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